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Abstract. Measurements of the individual multiplicities of , n~ and n produced in the deep-inelastic 
scattering of 27.5 GeV positrons on hydrogen are presented. The average charged pion multiplicity is 
the same as for neutral pions, up to 2 w 0.7, where z is the fraction of the energy transferred in the 
scattering process carried by the pion. This result (below z « 0.7) is consistent with isospin invariance. 
The total energy fraction associated with charged and neutral pions is 0.51 ± 0.01(stat.) ± 0.08 (syst.) and 
0.26 ± 0.01 (stat.) ± 0.04(syst.), respectively. For fixed z, the measured multiplicities depend on both the 
negative squared four momentum transfer and the Bjorken variable x. The observed dependence on 
agrees qualitatively with the expected behaviour based on NLO-QCD evolution, while the dependence on x 
is consistent with that of previous data after corrections have been made for the expected Q^-dependence. 
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1 Introduction 



The semi-inclusive production of pseudoscalar mesons in 
Deep-Inelastic Scattering (DIS) is a good tool to test the 
quark-parton model and QCD. A schematic diagram and 
the relevant variables for the process are shown in Fig. |l]. 
In the proton rest frame, the energy of the exchanged vir- 
tual photon ^* is v = E — E' {E and E' being the ener- 
gies of the incident and scattered positrons respectively), 
while its squared four-momentum is —Q^. The quantity 
X = /2Mv, where M is the proton mass, is the frac- 
tion of the light-cone momentum of the nucleon carried 
by the struck quark. The parton distribution function qf 
describes the momentum distribution of quarks in the nu- 
cleon, while the fragmentation function is a measure 
of the probability that a quark of flavour / fragments into 
a pion of energy Ej^ = zv. The quantity da/ is the cross- 
section for the absorption of the virtual photon by the 
struck quark. 




The quantity of interest in this paper is the pion differ- 
ential multiplicity, or the number (iV^) of pions produced 
in DIS, normalised to the total number (Nbis) of inclu- 
sive DIS events (e -l-p — > e' + X). In the QCD improved 
quark-parton model, it is given by the expression: 
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Y.e) f^dx qf{x,Q^)DJ{z,Q^ 



NoisiQ^) dz 



Ee/ Jodx (?/(a;,Q2) 



(1) 



Fig. 1. Semi-inclusive pion electroproduction diagram. 



where the sum is over quarks and antiquarks of flavour 
/, and e/ is the quark charge in units of the elementary 
charge. Perturbative QCD calculations in leading j^,^ and 
next-to-leading order |||,§,||j6[0,§,§ suggest a significant 
Q^-dependence of the fragmentation process. These QCD 
expectations have been verified by experimental results of 
hadron production in e'''e~ collisions (see references to ex- 
periments in R,pLpj), and in lepton-nucleon scattering at 
large v and ||ig,y • Data with pions 
identified in the final state are also available | |l8| , p^ , pO[j2 l| , p^ . 
This paper presents measurements at lower Q of multi- 
plicities for both charged and neutral pions as a function 
of the variable z in the range 0.1 < z < 0.9. This study 
of the Q^-dependence above 1 GeV^ provides new and 
precise information on the scaling violation of the pion 
fragmentation process at relatively small values of W (3.1 
GeV < W < 6.6 GeV) [||; 1^ is the invariant mass of 
the virtual-photon + proton system. 

Evidence for an additional x or W dependence of the 
multiplicities has been seen by previous experiments [^0|jO,|l^ . 
This has been ignored in the formalism leading to Eq.lU 
The multiplicities measured at HERMES were also stud- 
ied as a function of x to determine whether they show a 
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similar behaviour. Also, charged and neutral pion multi- 
plicities are compared as a test of isospin invariance. 

2 Experiment 

The measurements described here were performed with 
the HERMES spectrometer Q using the 27.5 GeV po- 
sitron beam stored in the HERA ring at DESY. The spec- 
trometer consists of two identical halves above and below 
the positron and proton beam pipes. The beam with a 
typical current in the range between 10 and 35 mA was 
incident on a hydrogen internal gas target [ p5| . The data 
used in this analysis were collected during the 1996 and 
1997 HERA beam periods. The target was operated in 
both unpolarised and longitudinally polarised configura- 
tions, with typical areal densities of 8x10^^ and 7x10^'^ 
atoms/cm^, respectively. The data from the polarised tar- 
get were analysed by averaging over the two spin orien- 
tations, producing results that are consistent with those 
from the unpolarised target. 

The scattered positrons and any resulting hadrons were 
detected simultaneously by the HERMES spectrometer. 
The geometrical acceptance of the spectrometer is ±(40 - 
140) mrad in the vertical direction and ±170 mrad in the 
horizontal direction. 

The identification of the scattered positron was accom- 
plished using a gas threshold Cerenkov counter, a transi- 
tion radiation detector, a scintillator hodoscope preceded 
by two radiation lengths of lead (prcshower counter), and 
an electromagnetic calorimeter. This system provided po- 
sitron identification with an average efficiency of 98% and 
a hadron contamination of less than 1%. Events were se- 
lected by imposing the kinematic restrictions > 1 GeV^ 
and y < 0.85, where y — v/E \s the virtual photon frac- 
tional energy. 

The Cerenkov counter was filled with a mixture of 70% 
nitrogen and 30% perfluorobutane (C4F10), providing a 
momentum threshold of 3.8, 13.6 and 25.8 GeV for pions, 
kaons, and protons respectively. The momentum was re- 
stricted to the range between 4.5 GeV and 13 GeV in the 
data analysis to eliminate possible kaon or proton con- 
tamination of the charged pion sample. The contribution 
from the inefficiency of the Cerenkov detector to the sys- 
tematic uncertainty on the charged pion multiplicities was 
evaluated to be 3.2% [|6). 

The electromagnetic calorimeter is composed of 840 
FlOl lead glass blocks |^^, and provided neutral pion 
identification by the detection of the two neutral clusters 
from two-photon tt*^ decay. Each of the two clusters was re- 
quired to have an energy > 1.4 GeV. The measurement 
of both the energies [E^i and E^2) and the relative angle 
{0.y^) of the two photons allowed the reconstruction of the 

invariant mass m~f^ — \J 4£'^ii?^2 sin^(6'-,-y/2). A typical 
measured spectrum of m^-y is shown in Fig. ^ The dis- 
tribution is centered at m^^ = 0.1348 ± 0.0008 GeV. The 
good tt" mass resolution of about 0.012 GeV allowed a safe 
background subtraction. The background was evaluated in 
each kinematic bin by fitting the invariant mass spectrum 



with a Gaussian plus a polynomial that reproduces well 
the shape of the background due to uncorrelated photons. 
The relevant contribution to the systematic uncertainty 
is less than 2%. It was evaluated by repeating the fitting 
procedure for different ranges and with polynomials of dif- 
ferent order . The number of tt*' detected was obtained 
by integrating the peak corrected for background, over the 
range ±2.5cr around the centroid of the Gaussian. 
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Fig. 2. Two-photon invariant mass spectrum. The solid line 
is a fit with a Gaussian plus a polynomial. The dashed line 
represents the background only. 



In order to exclude effects from nucleon resonances as 
well as kinematic regions with inadequate geometrical ac- 
ceptance, the additional requirement W"^ > 10 GeV^ was 
imposed for this analysis. This also helps to select hadrons 
originating from fragmentation of the struck quark, by 
excluding the kinematic region in which current and tar- 
get fragments are not well separated. Fragments from the 
target remnant are already strongly reduced due to the 
forward angle acceptance of the HERMES spectrometer. 

After background subtraction and data quality plus 
kinematic cuts, a total of 4. 2x10^ (1.3x10^) semi-inclusive 
charged (neutral) pions were considered in the analysis. 

The HERMES Monte Carlo program (HMC) was used 
to evaluate the detection probability for pions produced 
in DIS events. HMC is based on the LEPTO event gener- 
ator |2^ , the LUND fragmentation model ||2^ , and on the 
GEANT3 code for the simulation of the detector response 
[0. The LUND parameters were adjusted to fit vari- 
ous kinematic distributions from HERMES semi-inclusive 
data |^l[|. As an example. Fig. ^ shows comparisons be- 
tween the measured energy spectra of charged and neutral 
pions and the relevant HMC simulations. The agreement 
in the shape of the spectra is reasonable. The energy range 
of the charged pions is smaller than for neutral pions due 
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Fig. 3. Comparison between tiie measured and simulated 
pion energy spectra. Filled circles (triangles) are the neutral 
(charged) pion measurements. The histograms are the simu- 
lated spectra for neutral (dashed line) and charged (full line) 
pious. All spectra have been normalized to unit area. 



to the stronger kinematic restriction imposed on the data 
to ensure good charged pion identification. 
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Fig. 4. Detection probability as a function of z for charged 
(triangles) and neutral (circles) pions when the DIS positron 
is detected. 



This quantity was computed using HMC to account for ef- 
fects from pion losses due to the finite angular acceptance 
of the spectrometer, from detector inefficiencies, and from 
fiducial cuts imposed in the data analysis. A probability 
up to ^^40% is obtained for charged pions in the interme- 
diate and high z-region. At low z the detection probability 
for charged pions is further reduced due to the restriction 
on the pion momentum range. The probability for neutral 
pions is lower due to the need to detect both of the decay 
photons. Since the angle of these two photons is correlated 
with the tt" energy, the detection probability strongly in- 
creases with z, reaching ~18% in the highest z-bin. The 
positron detection probability does not affect the multi- 
plicity defined in Eq. |l]. Smearing by instrumental resolu- 
tion of the positron kinematics in the various v-hms was 
included in the Monte Carlo. This effect also leads to mi- 
gration of events in z. Corrections due to smearing vary 
from ~2% in the lowest z-bin to 30-35% in the highest z- 
bin. The systematic uncertainty on the multiplicities due 
to the detection probability (including smearing) is 7.5% 
(4.5%) for neutral (charged) pions. It was evaluated with 
HMC by drastically restricting the geometrical cuts or the 
allowed ranges of the relevant kinematic variables. The 
multiplicities plotted in this paper were corrected for the 
effects described above. 

Radiative corrections for internal bremsstrahlung pro- 
cesses were applied. These corrections range from ^3% in 
the lowest z-bin to ~15% in the highest z-bin, and result in 
a systematic uncertainty on the multiplicities of less than 
1%. Radiative corrections for external bremsstrahlung pro- 
cesses in the target were found to be negligible due to the 
small thickness of the hydrogen gas target |^ . 

It has been determined from Monte Carlo studies that 
diffractive production of pions is negligible in the HER- 
MES kinematic range, and so no correction was made for 
this contribution. 



3 Results 

Under the assumption of isospin invariance, the quark- 
parton model predicts that the multiplicity for neutral 
pions is equal to the average of those for positive and neg- 
ative pions. Specifically, the fragmentation function DJ 
is assumed equal to the average of the two charged pion 
fragmentation functions DJ and DJ , because the quark 
content of the vr*' is the same as the average of 7r+ and 7r~ 
The multiplicities and ^[^^i^ + ^]/2, 

are plotted as a function of z in Fig. B(a). Numerical val- 
ues are given in Tab. |^, where the individual multiplicities 
for 7r+ and 7r~ are also listed, along with the ranges and 
average values of the relevant kinematic variables. Both 
data sets in Fig. ||(a) show a strong decrease with the z 
variable that can be parameterized with the expression 



The detection probability of the pions produced in de- 
tected DIS events is shown as a function of z in Fig. ^. 



1 dN'' 



Ndis dz 



Nz"{l-zf. 
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The dashed curve shown in Fig. ||(a) is a indepen- 
dent fit to the tt" data using the above expression {N ^ 
0.335, a = -1.371, /? = 1.167). Also shown in Fig. |(a) 
is another Q^-independent parameterization using the in- 
dependent fragmentation model. The coefficients of this 
parameterization were obtained in Ref. by tuning on 
old data at Q^.^alues close to that of HERMES (see Q 
for references to these data). This parameterization re- 
produces the behaviour of the present data fairly well, 
apart from the high z-region. As expected from isospin 
invariance, the agreement between the data for neutral 
and charged pions is excellent, at least up to z ^ 0.70 as 
shown by the ratio in Fig. |^(b), which is consistent with 
unity below z ~ 0.70. The verification of isospin invariance 
agrees with the conclusion from earlier studies by BEBC, 
which used neutrino scattering 0,0, and by EMC |Q 
which, however, compared charged hadrons with neutral 
pions. The current data show that at higher z, the multi- 
plicity is larger for charged pions than for neutral pions. 
This difference suggests a possible contribution, via the 
radiative tail and instrumental resolution, from exclusive 
processes (e.g. 7* + p ^ 7r+ + Z\") where resonances af- 
fect each isospin channel differently. A recent calculation 
of hard exclusive electroproduction indeed predicts 
an enhancement of charged over neutral pion production 
of roughly an order of magnitude. Higher twist processes 
could also play a role at high z . The excess of charged 
over neutral pions at high z is being investigated further. 

Comparison of results obtained with unpolarised and 
polarised targets and of those with the top half and the 
bottom half of the spectrometer resulted in consistency 
within 3.5% and 2.5% respectively. It should be noted 
that the data sets for charged and neutral pions were 
obtained with very different event reconstruction proce- 
dures, detection efficiencies and background conditions. 
Hence the contributions to the systematic uncertainty are 
quite different for the two cases. The overall systematic 
uncertainty was estimated to be less than 9% (7%) for 
neutral (charged) pion electroproduction. The systematic 
uncertainty on the ratio plotted in Fig. ||(b) is 6%. 

The fraction of the energy v of the virtual photon 
transfered to pions 





.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



1 dN" 
' Ndis dz 



■dz 



(3) 



is 0.26 ± 0.01(stat.) ± 0.04(syst.) and 0.51 ± 0.01(stat.) ± 
0.08(syst.), for the neutral and charged cases, respectively. 
The corresponding number for neutral pions measured by 
EMC [H is 0.27±0.02±0.05. The current results indicate 
that the fraction of total energy carried by hadrons heav- 
ier than pions (mainly i^T"*" , p, p) is only ^ 23%. 
The integrals of Eq. ^ were evaluated by adding esti- 
mates of the contributions from the unmeasured z-regions 
to those from the measured z-region, the latter contribu- 
tions being 0.19±0.01(stat.)±0.017(syst.) for neutral and 
0.28 ± 0.01(stat.) ± 0.02(syst.) for charged pions. This ex- 
trapolation was based on the z-dependence given by the 
fit of ref. (solid line in Fig. H(a)). The total systematic 
uncertainty includes an estimate of the error in the extrap- 



Fig. 5. (a) Neutral (circles) and average charged (triangles) 
pion multiplicities. The error bars show the statistical uncer- 
tainty only. The systematic uncertainty for the charged (neu- 
tral) pions is 7% (9%). The solid line is a parameterization 
using the independent fragmentation model The dashed 
line is a fit to the present neutral pion data using the parame- 
terization given in the text. The charged pion data have been 
shifted slightly in 2 to make them visible, (b) Ratio of neutral 
to average charged pion multiplicities. The systematic uncer- 
tainty on the ratio (not included in the error bar) is 6%. 



6 



olation evaluated by comparing these results with those 
obtained using the fit of Eq. g. This contribution, 14% 
(13%) for charged (neutral) pions, dominates the overall 
systematic uncertainty. It is assumed in estimating this er- 
ror that the models used are a reasonable representation 
of the data, even at small z. The possibility of a radically 
different behaviour of the multiplicities at small z is not 
accounted for in the systematic uncertainty. 

In Fig. |(a) the HERMES results for the multiplicity of 
neutral pions as a function of z are compared with previ- 
ous results from EMC || and from SLAC Fig. |(b) 
compares HERMES results for the average charged pion 
multiplicity with fragmentation function data from EMC 
The HERMES results for both neutral and charged 
pions are systematically higher than those from EMC. 
This difference can be explained by the different range 
covered by the two experiments: (Q^) — 2.5 GeV^ for 
HERMES and (Q^) = 25 GeV^ for EMC. The range for 
the SLAC data is 1.8-8.5 GeV^. Similar QMependent be- 
haviour has been seen for hadrons in ref. |1^. In Figs. 0(a) 
and 0(b) the HERMES data have been evolved to the 
mean of the EMC data using an NLO model for the 
evolution of the fragmentation functions The agree- 
ment between the evolved HERMES data and the EMC 
data is much improved, especially for tt^, demonstrating 
the need for QCD corrections. The comparison of mul- 
tiplicities and fragmentation functions, and perhaps to a 
better extent the application to multiplicities of the model 
for the evolution of fragmentation functions, can be jus- 
tified if isospin symmetry is assumed (I? J = DJ; where 
TT represents the sum of positive and negative pions) and 
the strange quark contribution is neglected. In this case, 
the multiplicities are equivalent to the fragmentation func- 
tions. The effect of ignoring the strange quark can be in- 
ferred from a study in Ref. where the effect of ne- 
glecting all sea quarks was estimated to be less than 10% 
(between 20% and 40%) for favoured (disfavoured) frag- 
mentation functions. The dominance of it-quarks and the 
fact that the strange quark content is expected to be less 
than the light sea leads to the conclusion that the effect 
of the strange quarks should be well below 10%. On the 
other hand, the fact that the HERMES multiplicities are 
larger than the EMC fragmentation functions at low z, 
even after corrections (see Fig. 0(b)) could indicate a 
failure of the above approximation, which is more likely 
at low z. 

The charged pion multiplicities for HERMES are plot- 
ted in Fig. H as a function of x in four bins of z together 
with data on charged hadron multiplicities from EMC^ 
||ll|. AU data have been evolved to 2.5 GeV^, the 

average of HERMES, using the model described in 
the previous paragraph. The data are plotted at the mea- 
sured a;- value. The difference in the absolute value of the 
multiplicities is simply related to the fact that the HER- 
MES data are for pions only while the EMC data are for 
hadrons. A significant x-dependence is seen, which gets 
stronger as z increases. The mean per bin for HER- 
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Fig. 6^a) 7r° multiplicity from HERMES, EMC [gj and 
SLAC UM. (b) Average charged pion multiplicit y fr om HER- 
MES compared to EMC fragmentation functions [£0|. Only the 
statistical uncertainties are shown. The systematic uncertain- 
ties for neutral (charged) pions are 9% (7%) for HERMES, 
< 15% for SLAC and < 13% (< 10%) for EMC. 



^ The EMC pion data are not available in the required kine- 
matic binning. 
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Fig. 7. Multiplicities for (a) neutral pious, aud (b) the aver- 
age of charged pions (uote the EMC data for charged pions are 
for fragmeutation fuuctious). The HERMES results have beeu 
evolved to =25 GeV^ usiug a NLO QCD model §. Only 
statistical uncertainties are shown. The systematic uncertain- 
ties for neutral (charged) pions are 9% (7%) for HERMES and 
< 13% (< 10%) for EMC. 
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Fig. 8. Charged pion multiplicities (- 
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from HERMES (filled symbols) as a function of x in four dif- 
ferent 0-bins, compared to charged hadron multiplicities from 
EMC (open symbols). All data have been evolved to 2.5 GeV'^ 
and are plotted at the measured a;-value. 



MES varies only between 2.1 and 2.6 GeV'^. The observed 
x-dependence is therefore not generated by the correc- 
tion. It is striking that the slopes in the data from both 
experiments are consistent even though they were mea- 
sured at very different kinematics. The a;-dependence is 
not expected to affect the evolution of the data in- 
tegrated over the measured x range. This is bourne out 
by the good agreement after NLO-QCD corrections of the 
two data sets for 7r° multiphcitics in Fig. |7|(a). As stated 
earlier, the poorer agreement in Fig. ^(b) for charged pions 
could be due to the fact that multiphcitics are compared 
to fragmentation functions in this case. 

In order to better illustrate scaling violations in the 
fragmentation process, the Q^-dependence of the data at 
fixed z was studied. For this analysis, only the data in the 
range 0.03 < a; < 0.3 were considered, to allow for a use- 
ful range of Q^, which is quite limited by the HERMES 
acceptance for certain values of x. The total (neutral plus 
charged) pion multiplicity is plotted for 4 different z-bins 
in Fig. |[ The data show a clear Q^-dependence, espe- 
cially in the high z-bins. The variation of the data with 

(i.e. the slope of the curves) is mostly in agreement 
with the Q^-evolution predicted by the NLO QCD models 
[011:0 • It is worth mentioning that the above calculations 
are based on three analyses of fragmentation functions ex- 
tracted from e'^e~ data taken at center of mass energies 
of 29 GeV § and of 29 and 91 GeV |,|, i.e. at much 
higher energies than that of the data shown in Fig. ||. The 
sophistication of the parameterization increases for each 
successive version of the model, with the solid curve rep- 
resenting the most recent analysis. The absolute values of 
the present measurements are for the most part in rea- 
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Fig. 9. Total (neutral plus charged) pion multiplicity 
dQ^dz I "^^dcp-^ ^ function of for various z-bins. The sys- 
tematic uncertainty on the data is 8.5%. The three curves 
shown are NLO QCD calculations of fragmentation functions 
found in refs. (dashed line), D (dotted Hne), and §| (solid 
line) . 



sonable agreement with the QCD predictions, considering 
that each theoretical calculation has an uncertainty of up 
to 15%, the systematic uncertainty on the data is 8.5%, 
and that data on multiplicities are compared to a param- 
eterization of fragmentation functions. 



4 Conclusion 

Charged and neutral pion multiplicities in semi-inclusive 
deep-inelastic scattering at 27.5 GeV have been measured 
by the HERMES collaboration. The multiplicities are con- 
sistent with isospin invariance below z ~ 0.70, and show 
that about 3/4 of the energy transferred in the scattering 
is carried by pions. These measurements provide data at 
lower with improved statistical and systematic accura- 
cies compared to earlier measurements. The agreement of 
the current results with previous data is improved, espe- 
cially for tt", when an NLO Q^-evolution of the fragmen- 
tation process is taken into account. The Q^-behaviour 
of the present data resembles that of the NLO QCD cal- 
culations. An observed x-dependence of the multiplicities 
is similar to previous results on hadron production from 
EMC when the data are evolved to the same . 



5 Acknowledgments 

We thank M. Greco, G. Kramer, and B. Potter for use- 
ful discussions on fragmentation functions. We gratefully 
acknowledge the DESY management for its support and 



the staffs at DESY and the collaborating institutions for 
their significant effort. This work was supported by the 
Fund for Scientific Research-Flanders (FWO) of Belgium; 
the Natural Sciences and Engineering Research Council 
of Canada; the INTAS, HCM and TMR contributions 
from the European Community; the German Bundcsmin- 
isterium fiir Bildung, Wissenschaft, Forschung und Tech- 
nologic (BMBF) , the Deutscher Akademischer Austausch- 
dienst (DAAD); the Italian Istituto Nazionale di Fisica 
Nucleare (INFN); Monbusho, JSPS, and Toray Science 
Foundation of Japan; the Dutch Stichting voor Funda- 
menteel Onderzoek der Materie (FOM); the UK Particle 
Physics and Astronomy Research Council; and the US De- 
partment of Energy and National Science Foundation. 



References 

1. J.F. Owens, Phys. Lett. B76 (1978) 85. 

2. T. Uematsu, Phys. Lett. B79 (1978) 97. 

3. R. Baier and K. Fey, Z. Phys. C2 (1979) 339. 

4. N. Sakai, Phys. Lett. B85 (1979) 67. 

5. P. Chiappetta et al, Nucl. Phys. B412 (1994) 3. 

6. M. Greco, S. Rolh and A. Vicini, Z. Phys. C65 (1995) 277. 

7. J. Binnewies, B.A. Kniehl, and G.Kramer, Z. Phys. C65 
(1995) 471. 

8. J. Binnewies, B.A. Kniehl, and G. Kramer, Phys. Rev. 
D52 (1995) 4947. 

9. B. Kniehl, G. Kramer and B. Potter, Nucl. Phys. B582 
(2000) 514. 

10. EMC Collaboration, J.J. Aubert et ai, Phys. Lett. B114 

(1982) 373. 

11. EMC Collaboration, J. Ashman et al, Z. Phys. C52 (1991) 
361. 

12. P. Allen et al., Nucl. Phys. B176 (1980) 333. 

13. ZEUS Collaboration, J. Breitweg et al., Eur. Phys. J. Cll 
(1999) 251. 

14. HI Collaboration, C. Adloff et al., Nucl. Phys. B504 
(1997) 3. 

15. E665 Collaboration, M.R. Adams et al, Z. Phys. C76 
(1997) 441. 

16. P. Allen et al, Nucl. Phys. B214 (1983) 369. 

17. W. Wittek et al., Z. Phys. C40 (1988) 231. 

18. EMC Collaboration, J.J. Aubert et al., Z. Phys. C18 

(1983) 189. 

19. T.P. McPharlin et al., Phys. Lett. B90 (1980) 479. 

20. EMC Collaboration, M. Arneodo et al., Nucl. Phys. B321 
(1989) 541. 

21. EMC Collaboration, M. Arneodo et al., Z. Phys. C31 
(1986) 1. 

22. G. Drews et al., Phys. Rev. Lett. 41 (1978) 1433. 

23. P. Di Nezza, PhD thesis, Universita di Perugia (1998); 
HERMES internal note 98-092. 

(www-hermes. desy.de / notes /pub / doc-public-index.html) 

24. HERMES Collaboration, K. Ackerstaff et al., Nucl. Instr. 
and Meth. A417 (1998) 230. 

25. F. Stock et al., Nucl. Instr. and Meth. A343 (1994) 334. 

26. P. Geiger, PhD thesis, MPI Heidelberg (1998); HERMES 
internal note 98-005. 

(www-hermes. desy.de / notes /pub / doc-public-index.html) 

27. H. Avakian et al., Nucl. Instr. and Meth. A417 (1998) 69. 



28. G. Ingelman, A. Edin, and J. Rathsman, Comp. Phys. 
Comm. 101 (1997) 108. 

29. T. Sjostrand, Comp. Phys. Comm. 82 (1994) 74. 

30. GEANT Detector Description and Simulation Tool, CERN 
Program Library, Long Write-up W5014, 1994. 

31. H. Tallini, PhD thesis, University of Liverpool; HERMES 
internal note 98-024. 

(www-hermes. desy.de / notes /pub / doc-public-index.html) 

32. L Akushevich, N. Shumeiko, and A. Soroko, Eur. Phys. J. 
CIO (1999) 681. 

33. R.D. Field and R.P. Feynman, Phys. Rev. D15 (1977) 
2590. The data used in the fit corresponding to the solid 
line in Fig. ^ are shown in Fig. 6 of this reference. 

34. L.L. Frankfurt, M.V. Polyakov, M. Strikman, and M. Van- 
derhaeghen, Phys. Rev. Lett. 84 (2000) 2589. 

35. H. Avakian, Proceedings of the 8th International Work- 
shop on Deep Inelastic Scattering and QCD (DIS 2000), 
Liverpool, England, Apr 25 - 30, 2000. 



10 



z 


1 


1 iV""'*' 


1 i\ri-+ 


1 Tvr^ 




2 






0.142 


3.39±0.21 








0.223 


1 c\ci 1 c\ no 

1.9d±0.08 


o "I o_Ln "1 "1 

2.13±0.11 


z.zzzbU.UT 


^ r\ A _Ln no 

2.04±0.08 




1.49±0.05 


1.45±0.04 


1.63±0.02 


1.28±0.02 


0.324 


1.075±0.038 


1.0o±0.03 


1.22±0.0i 


0.897±0.023 


0.375 


0.774±0.028 


0.760±0.022 


0.906±0.017 


0.614±0.014 


0.424 


(J.5d3±U.Uz1 


U.55o±0.U2U 


ri C7n_i_rt r\"i cr 


n /I /io_i_n m o 

0.442±U.UlcJ 




n 491 -1-0 m fi 




n c;Qo-i_n ni i 
u.oozrnu.uii 




0.523 


0.342±0.014 


0.325±0.015 


0.411±0.012 


0.239±0.009 


0.574 


0.244±0.010 


0.246±0.011 


0.304±0.009 


0.188±0.007 
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0.192±0.010 


0.196±0.010 


0.250±0.009 


0.143±0.007 
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0.155±0.008 


0.161±0.009 


0.198±0.008 
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0.724 


0.115±0.007 
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0.164±0.007 


0.108±0.007 
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Table 1. Measured 7r°, {n'^ + n )/2, tv'^ , and tt , multiplicities for various 2;-bins. The average kinematic variables for this 
experiment arc: Etea^ = 27.5 GeV, ^ = 7.48 GeV, (Q^) = 2.5 GeV^ {W'^) = 28.6 GeV^ (u) = 16.1 GeV, (a;) = 0.082, while 
the ranges are [1, - 15] GcV^ W^= [10, ~ 44] GeV^ [~ 5.4, 23.4] GeV, and x= [0.03, ~ 0.6]. The quoted uncertainties 
are statistical. The systematic uncertainty for the charged (neutral) pions is 7% (9%). 



